SCIPP 94/28 
October 1994 



Cold + Hot Dark Matter Cosmology 
with m{i>ij) w m{i>r) ~ 2.4 eV 

'^ Joel R. Primack 

Qv^ Santa Cruz Institute for Particle Physics 

^""i University of California, Santa Cruz, CA 95064 

^ Jon Holtzman 

^ Lowell Observatory, Mars Hill Road, Flagstaff, AZ 86100 

•r^ Anatoly Klypin 

Department of Astronomy, New Mexico State University 

O Las Cruces, NM 88001 

(N 

O David 0. Caldwell 

Department of Physics, University of California 

^ Santa Barbara, CA 93106-9530 

CI, 
I 

O 



Abstract 



Cold + Hot Dark Matter (CHDM) Q = 1 cosmological models require a total 
neutrino mass ~ 5 eV. Because recent data support the z/^ — ?■ z/^ oscillation 
explanation of the cosmic ray z/^ deficit, which requires that ?n(z/^) ~ m(z/^), 
this suggests that ?n(z/^) ~ m(z/^) ~ 2.4 eV. The linear calculations and N- 
body simulation reported here indicate that an i7 = 1 CHDM model with 
two 2.4 eV neutrinos (designated Cz/^DM) agrees remarkably well with all 
available observations, but only if the Hubble parameter h ~ 0.5. We also 
show that even one 2.4 eV neutrino raises serious difficulties for low-i7 fiat 
CDM models. 
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Predictions of a Cold + Hot Dark Matter (CHDM) cosmological model with a single 
massive neutrino and ft^ = 0.3 (corresponding to m^ = 94/i^Oj, = 7 eV for Hubble parameter 
h = i7o/(100 km s~^ Mpc~ ) = 0.5) have been shown [1-3] to agree well with observations, 
with the possible exception that galaxies may form too late to account for the observations 
of quasars and damped Lya systems [4] at high redshifts ^; ^ 3 [5,6]. The latter observations 
can be accommodated [7] if the assumed u mass in CHDM is lowered from ~ 7 eV to ~ 5 
eV. Lowering the u mass in CHDM may also give a better account of the Void Probability 
Function [8] and of the properties of galaxy groups [3,9]; but with one ~ 5 eV u CHDM 
probably overproduces clusters, as we show below. 

Current experimental data suggest that the net u mass is shared among two species of 
neutrinos. Here we consider the consequences of such u masses for the formation of galaxies 
and large scale structure for cosmological models which are spatially flat and in which most 
of the DM is cold. We show that ft^ = 0.2 CHDM with the mass evenly shared between 
two u species — Ci/^DM — agrees better with observations than the one-neutrino version, 
better indeed than any other variant of CDM that we have considered. We first summarize 
the implications of the latest experimental hints regarding u masses, and then compare our 
calculations of cosmological predictions of various models with observations. 

It has been pointed out [10,11] that if the solar Ue and atmospheric u^ deficits arise 
from the existence of u masses, there are only two viable patterns for those masses: (A) 
the three active neutrinos are approximately degenerate; or (B) the nearly degenerate u^ 
and Ut constitute the hot dark matter, and the Ue and a sterile neutrino Us are lighter and 
are also nearly degenerate [12]. If one also takes into account the need for about 5 eV 
of neutrino mass for CHDM cosmological models to be viable, pattern (A) corresponds to 
m^^ sa m^^ sa m^^ sa 1.6 eV, while pattern (B) requires m^^ sa m^^ sa 2.4 eV. 

A u mass explanation of the solar Ue deficit, which is now fairly convincing, implies 
i^e -^ ^11 or Ue -^ ^s with i^ml- = |?7i(i/e)^ — ?7i(i/j)^| sa 10~^ eV^ between either pair of 
particles. Similarly, evidence for a u mass explanation of the deficit of i/^'s relative to 
i/e's in atmospheric secondary cosmic rays is also increasing, with compatible results from 
three experiments [13], and especially new information from Kamiokande [14]. The latter 
includes accelerator confirmation of the ability to separate Ue and u^ events, as well as an 
independent higher energy data set giving not only a u^/ue ratio agreeing with the lower 
energy data, but also a zenith-angle (hence source-to-detector) dependence compatible with 
i^fi -^ ^e or Uf^ — > Ut oscillations with '^m'l- sa 10~^ eV^. However, almost the entire region of 
A?7i^g — sin^ 2^^e allowed by the Kamiokande data is excluded by data from the Bugey and 
Krasnoyarsk reactor u oscillation experiments. Moreover, the absolute calculated Ue and u^ 
fiuxes — backed by measurements of /i fiuxes — agree with Ue data but show a u^ deficit 
[15]. Thus u^ — > Ut oscillations are favored as an explanation of the atmospheric u^ deficit. 

That the u mass pattern (B) might be correct is indicated by early results of the LSND 
experiment [16], the initial run of which showed an excess of about 8 beam-on events of a 
type which could be interpreted as i/^ — > i/g, whereas a background of ^ 1 event mimicking 
a Ue was expected. The LSND collaboration at this time is not claiming to have observed 
u oscillations, preferring to await results of their current run, which should give three times 
as much data. Nevertheless, the LSND positron energy distribution (if assumed to be from 
i^fi -^ 1^6 7 i^e -\- P ^ e'^ -\- n) appears [17] to be compatible with scheme (B) and not scheme 



(A), since the mass-squared difference required is A.m'^ ~ 6 eV^. This particuiar vaiue 
is not in conffict with the KARMEN [f8] or BNL E776 [f9] experiments, which are feast 
sensitive at that Am'^ where LSND is most sensitive. If the Ue mass is refativefy smaff (^ f 
eV, as indicated for Majorana u mass from neutrinoless double beta decay experiments), 
then the u^ mass is ~ 2.4 eV. This and the u^ — > Uj- explanation of the atmospheric u^ 
deficit then malces m[i'^) sa ?7i(i/^) sa 2.4 eV. It is this scenario for the hot dark; matter 
in a CHDM cosmology which we will show below gives predictions that are in remarkable 
agreement with astronomical observations. 

COBE observations [20] of fiuctuations in the microwave background radiation provide 
an upper limit (since they include possible tensor gravity wave as well as scalar density wave 
contributions) on the normalization of the spectrum of fiuctuations. When this normal- 
ization is used for the standard Cold Dark Matter (CDM) model [21] in a critical density 
(0 = 1) universe with a Zel'dovich primordial power spectrum {P{k) = Ak"'P with Up = 1) 
as predicted by simple infiationary models, this fits large-scale data but produces too much 
structure on smaller scales; for example, there are far too many clusters of galaxies [22] and 
small-scale velocities of galaxies are too large [3,9]. 

CDM is attractive because of its simplicity and the existence of well-motivated particle 
candidates (lightest superpartner particle and axion [23]) for the cold DM; moreover, CDM 
came remarkably close to predicting the COBE signal. So several variations have been tried 
to patch up the CDM model. Lowering the normalization (introducing "bias") or "tilting" 
the primordial spectrum (assuming Up sa 0.7) improves agreement with data on small scales 
at the cost of poorer agreement on large scales. The variants of CDM [24] that agree best 
with observations add either a cosmological constant (ACDM) or a little hot (neutrino) dark 
matter (CHDM). 

We report (quasi-)linear estimates for various observable quantities in Table I. All models 
in the Table are normalized to COBE except for the CDM model labelled "biased". Actually, 
our COBE normalization [Q = 17 /iK) is about 10% lower than the latest analyses [20] would 
suggest. We have chosen this normalization to allow for a little gravity wave contribution and 
tilt (consistent with the expectations from simple models of cosmic infiation), and we regard 
it as being both more realistic and more conservative than a higher normalization — more 
conservative, since the greatest problem for the CHDM models is enough early structure 
formation, a problem that worsens as the normalization of the spectrum is decreased. 

The first two lines of numbers give our estimates of a variety of observational quan- 
tities and the uncertainties in them, from large to small scales. The bulk velocity at 
r = 50 h~^ Mpc is derived from the latest POTENT analysis [25]; the uncertainty includes 
the error from the analysis but not cosmic variance. However, similar constraints come from 
other data on large scales such as power spectra that may be less affected by cosmic variance 
since they probe a larger volume of the universe. We have estimated the current number 
density of clusters (A^ciust) from comparison of data on the cluster temperature function from 
X-ray observations with hydrodynamic simulations [26] as well as from number counts of 
clusters [27]. All recent estimates of the cluster correlation function give fairly large values 
at 30 h~^ Mpc [28]; this also suggests that the zero crossing of the correlation function must 
exceed ~ 40 h~^ Mpc. The linear estimate of pairwise velocities (ay) is not an observed 
value, since pairwise velocities are strongly influenced by nonlinear evolution. However, from 
experience with N-body simulations for various models, we have found that the results from 



simulations are roughly a factor of 3 larger than the linear estimate. We have confirmed 
that factor by such a simulation for the Ci/^DM model (high resolution 800"^ PM mesh in a 
50 h~^ Mpc box with 256"^ cold and 2 x 256"^ hot particles), but all values given here are for 
linear calculations. The limit we give is therefore our estimate of the maximum linear value 
allowed by observations. The final column gives the observed density in cold hydrogen and 
helium gas at z = 3.0 — 3.5 from observations of damped Lyman a systems [4]. 

The next two lines present predictions from the CDM model, and illustrate its problems. 
The cluster correlation function at 30 h~^ Mpc is smaller than observations indicate regard- 
less of CDM normalization, refiecting the fact that the matter correlation function becomes 
negative beyond 36 h~^ Mpc. If CDM is normalized to erg = 0.7 (or equivalently to linear 
bias b = ag^ = 1.43), the cluster density problem is avoided, but small-scale velocities are 
still too large and bulk velocities on a scale of 50 h~^ Mpc are probably too low. 

The problem with CDM is that it has too much power on small scales relative to power 
at large scales. Since including hot DM reduces small scale power (because u free streaming 
causes perturbations to damp on smaller scales), including a u component improves the 
agreement with observations. The model studied by KHPR [1] with a u mass of 7 eV, 
corresponding to ft^ = 0.3 for h = 0.5, is a much better match to observations than CDM, 
but it has Ogas too small [6,7]. The small-scale velocities in this model are small enough 
[1] to agree with the old result cry[l h~^ Mpc) = 340 km s~^ from the CfAl survey [29]. 
However, this result is now known to be in error because of the accidental omission of 
the Virgo cluster [30]; moreover, the (Ty statistic is not very robust [31], since it is heavily 
infiuenced by the presence of (relatively rare) clusters. A direct comparison of galaxy groups 
in "observed" CDM and CHDM simulations with identically selected CfAl groups shows 
that CDM velocities are much too high, even with biasing, while the velocities in the ft^ = 0.3 
CHDM model are in reasonable agreement [3,9]. However, the fraction of galaxies in groups 
is slightly too high for ft^ = 0.3 CHDM, while it is significantly too low for CDM. Thus 
agreement is improved for a lower ft^. 

Lowering ft^ to 0.20 (lu) increases small-scale velocities but not enough to confiict with 
the data, and it raises Ogas enough for early object formation [7]. However, this model 
probably overproduces clusters. In order to avoid this, it could be normalized lower, given 
some gravity wave contribution to COBE, but this would result in too little early structure 
formation. 

Now consider Ci/^DM. All quantities are in good agreement with the astronomical data 
if the same ft^ = 0.2 is divided between two u species, as suggested by the data. The ratio 
of the power spectrum for Ci/^DM compared to that for CHDM with the same total u mass 
in one species is ~ 1 at large and small scales, but it has a dip of about 30% centered at 
~ 10 h~^ Mpc. The larger u free-streaming length, resulting from a u mass of 2.4 eV instead 
of twice that, lowers the abundance of clusters and agrees better with observations. 

It is remarkable that, with the experimentally suggested u masses, only cosmological 
models with h sa 0.5 match observations. Returning to Table I, note that for h = 0.7 — 
favored by many observers — CDM (CDMq.t) is an even worse fit to the data than for 
h = 0.5 because the larger h makes matter-dominance (oc fth^) occur earlier and thus moves 
the bend in the CDM spectrum to smaller scales, giving more intermediate and small scale 
power for a given large scale normalization. Adding two 2.4 eV neutrinos (Ci/^DMq.t) only 
slightly improves the situation, because this only gives ft^ ~ 0.1 for h = 0.7, so the spectrum 



is not modified very much. (Recall that for given m{u), Oj, scales as h~^ since critical density 
is oc h^.) The h = 0.7 model can match data better with a tilted spectrum of primordial 
fluctuations with Up = 0.8f (Cu^DMno.s) — but only without gravity waves. Typical cosmic 
inflation models with this much tilt would suppress scalar power by a factor of ~ 0.5 leading 
to very serious underproduction of clusters and of Ogas, although models have been proposed 
[32] with tilt but no gravity waves. Of course, with large /i, = 1 leads to too short a time 
since the big bang: to = f-ffo"^ = 6.52 Gy/h = 9.3 Gy for h = 0.7. 

A larger age is obtained for Oq < 1 which, to be consistent with inflation, requires 
a positive cosmological constant A. The maximum value allowed by the COBE data is 

Oa = ^/(3-ffo) ~ 0.78 [33], and the maximum allowed by quasar lensing statistics is Oa ~ 0.7 
[34]. For a flat [k = 0) universe with Oa = 0.7 and Oq = 0.3, h = 0.7 corresponds to to = 13.5 
Gy. ACDM with these parameters is a fairly good fit [35] to the data, although our linear 
calculations suggest that not enough clusters are produced and bulk velocities may be too 
low. However, this model becomes much worse if even one u of 2.4 eV is added, seriously 
underproducing clusters and Ogas- 

A similar situation occurs for = 1 Ci/^DM with h = 0.4, for which Oj, = 0.32 with 
two 2.4 eV neutrinos. Because the bend in the CDM spectrum moves to larger scales as h 
decreases, there is less small scale power for given large scale normalization; adding hot DM 
further decreases small scale power. We find that even with only one 2.4 eV i/, there is just 
not enough power to generate the observed number of clusters or high-redshift objects. 

Ever since the early 1980s there have been hints [36] that features on small and large 
scales may require a hybrid scenario in which there are two different kinds of dark matter. 
Preliminary studies of the CHDM scenario were carried out in 1984 [37], and it was first 
worked out in detail only in the last two years [1-3] with one massive u. We have shown 
here that the Ci/^DM model, with Hubble parameter h = 0.5 and both neutrinos having 
a mass of 2.4 eV as suggested by ongoing experiments, gives a remarkably good account 
of all presently available astronomical data. New data on CMB, large scale structure, and 
structure formation will severely test this highly predictive model. Experimental results 
expected soon will clarify whether indeed m[i'^) sa ?7i(i/^) sa 2.4 eV. Table I shows the 
implications of such u masses for a variety of popular CDM-type cosmological models. If 
even just the u^ has a mass of 2.4 eV, as suggested by preliminary results from the LSND 
experiment, fiat low-0 CDM models are disfavored. 
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TABLES 

TABLE L Comparison of models: COBE normalization (Qrms-ps = 17 /^K 
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ACDM/ACHDM models. 
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^Niy is the number of v species with mass. If Aj, > 1, each species has the same mass mj,. 
b(AM/M)™s for i?top-hat = 8/1-1 ^pc. 
■^Bulk velocity in top-hat sphere of radius 50h~^ Mpc. 

^^Number density of clusters A(> M) in units of W~'^ h'^yLpc~'^ above the mass M = W^^h~^MQ, 
calculated using Press-Schechter approximation with gaussian filter and Sc = 1.50. 
''The cluster-cluster correlation function amplitude at 30h~^ Mpc, computed using linear theory 
[24] and assuming a unit bias factor for the dynamical contribution. 
Zero crossing (^(r) = 0) of the correlation function in units of h~^ Mpc. 

^Linear estimate of pairwise velocity at r = lh~^ Mpc scale: a^ = 2Hq J dk P{k) (1 — sin kr)/kr. 
'^Mean density of collapsed baryons at 2: = 3 — 3.5 in units of 10-"^ of critical density, calculated using 
rigas = (rife/r^c) erfc(5c/'v2'7), with 5c = 1.4 [7], and a computed for mass 5 X W^^h~^MQ using 
gaussian smoothing and assuming all gas is neutral. Since some gas may be ionized or removed by 
star formation, i7gas for the various models should be at least as high as the observations. 



